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The strong CP problem was solved by Peccei Sz Quinn by introducing axions, which are a viable candidate for 
DM. Here the PQ approach is modified so to yield also Dark Energy (DE), which arises in fair proportions, 
without tuning any extra parameter. DM and DE arise from a single scalar field and, in the present ecpoch, 
are weakly coupled. Fluctuations have a fair evolution. The model is also fitted to WMAP release, using a 
MCMC technique, and performs as well as LCDM, coupled or uncoupled Dynamical DE. Best-fit cosmological 
parameters for different models are mostly within 2-cr level. The main peculiarity of the model is to favor high 
values of the Hubble parameter. 


1. Introduction 

Axions are likely to be the Dark Matter (DM) that 
cosmological data require. Axions arise from the so¬ 
lution of the strongCP problem proposed by Peccei 
& Quinn in 1977 ([j, PQ hereafter), who suggested 
that 9 parameter, in the effective lagrangian term 

Ce = ^9G-G ( 1 ) 

ZTT 

{as', strong coupling constant, G: gluon field ten¬ 
sor), causing GP violations in strong interactions, is 
a dynamical variable. Under suitable conditions 9 ap¬ 
proaches zero in our epoch, so that the term m is sm- 
pressed, while residual 0 oscillations yield DM [3,|3|- 

In the PQ scheme, 9 is the phase of a complex field 
$ = its evolution is set by the potential 

U(|$|) = A[|<I>P-F|qP , (2) 

whose U{1) invariance breaks when T < Fpg (the 

PQ energy scale, which shall be ~ lO^^GeV). Then 
(j) settles at the potential minimum, while 9 takes a 
random value, different in different horizons. When 
chiral symmetry also breaks down, during the quark- 
hadron transition, a further term 

Ui = [^(0(T)|TO|0(T))m,] (1 - cos9) (3) 

Q 

must be added to the lagrangian, because of instan- 
ton effects. At T ~ 0, the square bracket approaches 
{mTT, fn- TT-meson mass, decay constant). 

The dependence on the site of initial 9’s will 
causes later adjustments, as soon as the potential © 
switches on. When 9 is small, Vi oc 9"^ and the 9 field, 
undergoing harmonic oscillations, is DM. 

Here we wish to discuss recent work [4], where the 
NG potential (0 is replaced by a tracker potential 0. 
Then, instead of settling on a value Fpq , (j) continues 
to evolve over cosmological times, at any T. As in 
the PQ case, the potential involves a complex field $ 
and is U{1) invariant. Its phase is horizon-dependent 


but, just as in the PQ case, evolves only after chiral 
symmetry breaks, yielding a term similar to ®. 

At variance from the PQ case, however, the 9 evolu¬ 
tion starts and continues while also (j) is still evolving. 
This goes on until our epoch, when </> is expected to 
account for Dark Energy (DE), while, superimposed 
to such slow evolution, faster transversal 9 oscillations 
occur, accounting for DM. This scheme, where the 4) 
fields bears a dual role, with its modulus and phase, 
will be dubbed here dual axion model. 

The question is whether (j) and 9 dynamics interfere, 
still allowing 9 to take values so small to solve the CP 
problem and preserving harmonic oscillations allow¬ 
ing the prescribed amount of DM. Here we show that 
the reply is positive. As it can be expected, however, 
DM and DE are dynamically coupled, although this 
coupling weakens as we approach the present era. 

Here we also consider a generalization of the 
SUGRA potential as an example of tracker poten¬ 
tial. With an energy scale A ^ 10^°GeV, this model 
allows suitable values for today’s DM and DE densi¬ 
ties, while 9 is even smaller than for PQ. DM density 
fluctuations are also able to grow and to account for 
the observed large scale structure. We fit the expected 
angular GBR spectra to data, finding that this model 
substantially performs as well as ACDM. 


2. Lagrangian theory 

In the dual-axion model we start from the la¬ 
grangian C = - U(|4)|)} which can 

be rewritten in terms of (j) and 9, adding also the term 
breaking the U{1) symmetry, as follows: 

— U((/)) — m^(r, 0)(/)^(l — COS0)} . 

Here is the metric tensor. We shall assume that 
ds^ = Pfj^ijdx^dx'^ = af{dT‘^ — rjijdxidxj), so that a is 
the scale factor, r is the conformal time; greek (latin) 
indeces run from 0 to 3 (1 to 3); dots indicate dif¬ 
ferentiation in respect to r. The mass behavior for 
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T ~ Aqcd will be detailed in Section 3. The equa¬ 
tions of motion read 

9-\-2{a/a-\-(j)) 9-\-wfaf sinO — Q , (4) 

(j)-\-2{a/a) ^a^V{(f) = (fO"^, (5) 

while the energy densities = Pe,(t>-kin +Pe,<t>-pot and 
pressures= pe,<j>-kin-Pe,(i>-pot, under the condition 
0 <C 1, are obtainable from 

P0,kin = (f^9^ 120? , pe,poi = m^(T, , 

P(jy,kin — ^ /2a , P(j^,pot — ^ (0) • (6) 


3. Axion mass 

According to eq. 0 . the axion field begins to os¬ 
cillate when m(T, (/))a ~ 2( a/a -|- (/>/(/>). In the dual- 
axion model, just as for PQ, axions become massive 
when the chiral symmetry is broken by the formation 
of the qq condensate at T ~ ^qcd- Around such T, 
therefore, the axion mass grows rapidly. In the dual- 
axion model, however, a slower growth takes place also 
later, because of the evolution of (j). Then m(T, </) is 

mo((/) = mT,fT,/cj) = (0.0062/^) GeV . (7) 

Since (j) ~ rup today, the axion mass is now TOq ~ 
5 • I0“^^eV, while, according to at high T : 

m(r,0)~O.lmo((/)(AQCD/T)3® . (8) 

This expression must be interpolated with eq. 0 , to 
study the fluctuation onset for T ~ Aqcd- We solved 
the equations of section 2 by assuming 

m{T,(j)) = a < a, 

m{T, (p) = mo{4>) a > Gc 

with Gc = 2.16 • 10“^^. With the selected value, when 
T < 0.5 Aqcd, m{T,(j)) already approaches its low-T 
behavior mo{(j>). Let us finally outline that eqs. (TJ 
and © , as well as the above interpolation, show that 
m(T, (/)(/) is (/-independent, as is required to obtain 
the equation of motion 0. Fig.0 shows the (low-)^ 
dependence of mo{(j)). Notice the rebounce at z ~ 10, 
whose impact on halo formation could be critical |3|- 


4. The case of SUGRA potential 

When 9 performs many oscillations within a Hubble 
time, then (pe^kin) — {pe,pot) and (pg) ~ 0. By using 
eqs. 0 , 0 ,®, it is easy to see that 

771 777 

pe + iHpg = — pe , p0-l-3i7(p0-l-p0) =- pg . (9) 



log(a/a„) 

Figure 1: p variations cause a dependence of the effective 
axion mass on scale factor a, which is shown here. 

When m is given by Eq 0 , m/m = —(//(/ —3.8 T/T. 
At r ~ 0, instead, m/m ~ —<//(/. Here below, the 
indices g, cj, will be replaced by dm,de- Eqs. 0 show 
an energy exchange between DM and DE. The former 
eq. 0 can then be integrated, yielding pdm cx m/cff. 
This law holds also when T <C Aqcd, and then the 
usual behavior pdm ok a~^ is modified, becoming 

PDMa^/) — const. (10) 

Let us now assume that the potential reads 

1/(0) = (A“+^/(/“)exp(47r(/Vwp) (11) 

(no 9 dependence); in the radiative era, it will then be 

0“+2 = 5„A“+4aV" , (12) 

with Pa = a{a 2)^/4(q; -|- 6). This tracker solution 
holds until we approach the quark-hadron transition. 
Then, in Eq. 0, the DE-DM coupling term, 00^, 
exceeds a^V and we enter a different tracking regime. 

This is shown in detail in Fig. 0, obtained for mat¬ 
ter (baryon) density parameters 11^ = 0.3 {Dg = 
0.03) and Hubble constant h = 0.7 (in units of 100 
km/s/Mpc). In particular, Fig. 0l) shows the tran¬ 
sition between these tracking regimes. Fig. 0) then 
shows the low-z behavior (l-|-z = 1/a), since DE den¬ 
sity exceeds radiation (z ~ 100) and then overcomes 
baryons (z ~ 10) and DM (z ~ 3). Fig. 0) is a 
landscape picture for all components, down to a = 1. 

Notice the a dependence of pdm, occurring accord¬ 
ing to Eq. (irnil . In Fig. 0 the related behaviors of the 
density parameters fli {i = r, b, 9, 0, i.e. radiation, 
baryons, DM, DE) are shown. 

In general, once the density parameter Ude (at 
z = 0) is assigned, a model with dynamical (coupled or 
uncoupled) DE is not yet univocally determined. For 
instance, the potential CD depends on the parame¬ 
ters OL and A and one of them can still be arbitrarily 
fixed. Other potentials show similar features. 

In this model such arbitrariness no longer exists. 
Let us follow the behavior of pdm, backwards in time. 
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Figure 2: Density behaviors vs. scale factor. 

until the approximation 0 <C 1 no longer applies. This 
must approximately coincide with the time when 

2{a/a-\-^/(j)) ~ m{T, (f) a (13) 

and 9 enters the oscillation regime. At that time, ac¬ 
cording to Eq. m, which is marginally valid up to 
there, and taking 9=1, 

^ ^ </*(«)] (a) • (14) 

(p{a) a'^ 

The system made by eqs. 113 and d, owing to 
eq. m, allows to obtain (i) the scale factor Oh when 
fluctuations start and (ii) the scale A in the potential 
d- To do so, the present density of DM, Po.DM, 
(or Ddm) must be assigned. But, as we shall soon 
outline, the observational value of the density in the 
world forces Oh to lay about the quark-hadron tran¬ 
sition, while also A is substantially fixed. 

The plots in the previous section are for Udm = 
0.27; then A ~ 1.5 • 10^°GeV and an ~ 10“^^ are 
required by the system. But, when Ddm goes from 
0.2 to 0.4, logiQ(A/GeV) (almost) linearly runs in the 
narrow interval 10.05-10.39 , while an steadily lays at 
the eve of the quark-hadron transition. 



Figure 3: Density parameters fli vs. scale factor a. 

5. Evolution of inhomogeneities 

Besides of predicting fair ratios between the world 
components, a viable model should also allow the for¬ 
mation of structures in the world. 

The dual axion model belongs to the class of cou¬ 
pled DE models treated by Amendola @ , with a time- 
dependent coupling = Ijf). A (j)-MDE phase 
therefore exists, after matter-radiation equivalence, as 
the kinetic energy of DE is non-neglegible during the 
matter-dominated era. 

Fluctuation evolution is then obtained by solving 
the equations in 0, with the above C'((()). The behav¬ 
ior shown in Fig.Q) (left) is then found. 

Fig 0 (center and right) compare fluctuation evo¬ 
lutions the dual axion model (solid curves), with those 
in an analogous AGDM model (dot-dashed curves) 
and in a coupled DE model with constant coupling 
C = 0.25 VSttG ~ {C{4>)) (dashed curves). As shown 
by the plots, the overall growth, from recombination 
to now is similar in dual axion and AGDM models, 
being quite smaller than in DE models with constant 
coupling. The differences of dual axion from AGDM 
are: (i) objects form earlier and (ii) baryon fluctua¬ 
tions keep below DM fluctuations until very recently. 


6. Comparison with WMAP data 

We shall then test the dual axion model towards 
GMB data, together with other dynamical or coupled 
DE cosmologies, by using a parameter space of 7 to 
8 dimensions. We use a Markov Chain Monte Carlo 
(MCMC) approach (e.g. i), just as in the original 
analysis of WMAP first-year data 0, constraining 
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Figure 4: Time evolution of DM and baryon fluctuations. 


flat ACDM in respect to six parameters: 
h, n, the fluctuation amplitude A and the optical 
depth r. Notice that, with the naming convention 
used here, Vim = + VIum- 

Here, three classes of DE are considered: (i) 
SUGRA dynamical DE, requiring the introduction of 
the parameter A = log^Q(A/GeV), the energy scale 
in the potential (12). (ii) Constant coupling SUGRA 
DE, requiring a further parameter, the coupling (3 = 
C (3mp/167r)^/^. ( iii) Coupled models including the 
dual axion case. For this model the parameter /3 does 
not exist, being C = (f>~^ and also the scale A is con¬ 
strained by the requirement that VI^e lays in a fair 
range (also solving the strong CP problem), so that A 
and VIee are no longer independent parameters. The 
(iii) class of model, that we call models, is how¬ 
ever set keeping A as a free parameter. We aim to 
test whether data constrain A into the right region, 
turning a generic model into a dual-axion model. 

To use the MCMC, we need an algorithm provid¬ 
ing the Cl's. Here we use our optimized extension of 
CMBFAST able to inspect the cosmologies (i), 
(ii) and (iii). Then, the likelihood of each model is 
evaluated through the publicly available code by the 
WMAP team ^ and accompanying data [l^ . 

6.1. Results 

The basic results of our analysis are summarized in 
the Table |3 El and IHII For each model category we 
list the expectation values (x) of each parameter x 
and the associated variance Ux] we also list the values 
of the parameters of the best fitting models Xmax- 

Let us however remind that 4)~^ models, whose fit¬ 
ting results are reported in Table lHlI include the dual- 
axion model, but many other cases as well. Our ap¬ 
proach was meant to test whether CMB data carry 
information on A and how this information fits the A 
range turning a 4)~^ model into the dual-axion model. 

At variance from uncoupled and costant-coupling 
SUGRA models, WMAP data yield constraints on A 
for models (see Table IHTll and the 2-cr A-interval 
ranges from ^ 10 to ^ 3 • 10^°GeV, so including the 
dual-axion model. In Fig. 0 the Cf spectra for all 
best-fit models (apart of ACDM) are compared. At 


Table I SUGRA parameters for uncoupled DE 


X 

{x) 

O'x 

Xmax 


0.025 

0.001 

0.026 

VloMbf 

0.12 

0.02 

0.11 

h 

0.63 

0.06 

0.58 

T 

0.21 

0.07 

0.28 

ris 

1.04 

0.04 

1.08 

A 

0.97 

0.13 

1.11 

A 

3.0 

7.7 

13.7 


Table II SUGRA parameters in the presence of a 
constant DE-DM coupling j3. 


X 

(x) 

O'x 

Xmax 

Vlhh^ 

0.024 

0.001 

0.024 

VIdmIS 

0.11 

0.02 

0.12 

h 

0.74 

0.11 

0.57 

T 

0.18 

0.07 

0.17 

ris 

1.03 

0.04 

1.02 

A 

0.92 

0.14 

0.93 

A 

-0.5 

7.6 

8.3 

/3 

0.10 

0.07 

0.07 


Table III SUGRA parameters for & <p ^ model. 


X 

{x) 

Ox 

Xmax 


0.025 

0.001 

0.026 

VloMh? 

0.11 

0.02 

0.09 

h 

0.93 

0.05 

0.98 

r 

0.26 

0.04 

0.29 

Us 

1.23 

0.04 

1.23 

A 

1.17 

0.10 

1.20 

A 

4.8 

2.4 

5.7 


large I all models yield similar behaviors and this is 
why no model category prevails. Discrimination could 
be achieved by improving large angular scale observa¬ 
tion, especially for polarization, so to reduce the errors 
on small-Z harmonics. 


7. Discussion 

A first point worth outlining is that SUGRA (un¬ 
coupled) models, bearing precise advantages in re¬ 
spect to ACDM, are consistent with WMAP data. 
The ratio w = p/p at z = 0, for most these mod¬ 
els, fulfills the constraint w ^ —0.80. However, these 
models exhibit a fast w variation, and w becomes ^ 
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Figure 5: Cf spectra for the best fit SUGRA (solid line), 
constant coupling (dotted line), ())“^-coupling (dashed) 
and dual-axion (dot-dashed) models. 


-0.6, -0.7 at 2 ~ 1-2. In spite of this sharp decrease, 
however, there is no conflict with data. 

The results shown in the previous section were also 
considered in the presence of some priors. For un¬ 
coupled or constant-coupling SUGRA models, adding 
priors scarsely affects conclusions. ())“^-models, in¬ 
stead, are to be discussed separately. 

In the first two categories of models opacity (r) is 
pushed to values even greater than in ACDM (see 
also 0 ). This can be understood in two comple¬ 
mentary ways: (i) Dynamical DE models, in general, 
exhibit a stronger ISW effect, increasing Cf in the 
low-Z plateau (e.g. ^3)- To compensate this effect, 
in the fit of data, the spectral index rig is pushed to 
greater values. In turn, owing to the r-n degeneracy, 
this is compensated by increasing r. (ii) If the TE 
correlation is also considered, it is worth reminding 
that dynamical DE lowers the TE correlation at low 
I ^ 3 . To fit the same observed correlation level, a 
greater r is therefore favored. 

Greater t’s have an indirect impact also on 
whose best-fit value becomes greater, although con¬ 
sistent with ACDM within 1-cr. If we impose, as a 
prior, = 0.0214±0.0020 (BENS estimates [13), 
h shifts slightly below HST findings, still well within 
1-(T. We shall consider the effect of a prior also on h. 

The prior on D},h^ affects reionization and the spec¬ 
tral slope: r and rig are lowered to match WMAP’s 
findings. The prior on h, in the absence of coupling, 
favors low-A models, closer to ACDM. Quite in gen¬ 
eral, in fact, the sound horizon at decoupling is unaf¬ 
fected by the energy scale A, while the distance of the 
last-scattering-band is smaller for greater A’s. Then, 
as A increases, a lower h is favored to match the po¬ 
sition of the first peak. In the presence of coupling, 
there is a simultaneous effect on /3: greater /3’s yield 


a smaller sound horizon at recombination, so that the 
distribution on h is smoother. 

0“^ models exhibit rather different features. Pa¬ 
rameters are more strongly constrained in this case, 
as already outlined for the scale A. The main puz¬ 
zling feature of models is that large h is favored: 
the best-fit 2-a interval does not extend below 0.85. 

The problem is more severe for dual-axion A’s val¬ 
ues, whose preferred values are well within 2-a. This 
model naturally tends to displace the first Cf peak 
to greater I (smaller angular scales) as coupling does 
in any case does. The model, however, has no ex¬ 
tra coupling parameter and the intensity of coupling 
is controlled by the scale A. Increasing this scale re¬ 
quires a more effective compensation and still greater 
values of h are favored. 

This effect appears however related to the choice 
of SUGRA potentials, which is just meant to pro¬ 
vide a concrete framework for the dual axion model. 
To deepen its analysis, the contribution to the axion 
abundance due to the decay of possible topological 
defects should also be discussed. 

A previous analysis of WMAP limits on constant 
coupling models had been carried on by M Their 
analysis concerned potentials V fulfilling the relation 
dV/df) = BV^, with suitable B and N. They also 
assumed that r = 0.17. Our analysis deals with a 
different potential and allows more general parame¬ 
ter variations. The constraints on j3 we find are less 
severe. It must be however outlined that [3 > O.I- 
0.2 is forbidden by a non-linear analysis of structure 
formation Q- 

8. Conclusions 

The first evidences of DM date some 70 years ago, 
although only in the late Seventies limits on GMB 
anisotropies made evident that a non-baryonic com¬ 
ponent had to be dominant. DE could also be dated 
back to Einstein’s cosmological constant, although 
only SNIa data revived it, soon followed by data on 
GMB and deep galaxy samples. 

Axions have been a good candidate for DM since the 
late Seventies, although various studies, as well as the 
occurrence of the SN 1987a, strongly constrained the 
PQ scale around values 10^° < Fpg < lO^^GeV. Gon- 
tributions to DM from topological singularities (cos¬ 
mic string and walls) narrowed the constraints to Fpg. 
Full agreement on the relevance of such contributions 
has not yet been attained and, in this paper, they are 
still disregarded, while they could cause shifts in our 
quantitative predictions. This point must be therefore 
deepened in further work. 

The fact that DM and DE can both arise from scalar 
fields, just by changing the power of the held in effec¬ 
tive potentials, already stimulated the work of various 
authors. A potential like CD was considered in the 
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Figure 6: The onset of coherent axion oscillations, due to 
the increase of m{T, <j>), causes the behaviors of pe,pot, 
p9,kiri (a) and =p,p,e/p ^,0 (b) shown here. 

so-called spintessence model [^. According to the 
choice of parameters, $ was shown to behave either 
as DM or as DE. Padmanabhan & Choudhury 
instead, built a tachionic model where DM and DE 
arise from a single scalar field. 

Here, the scalar field <&, accounting for both DE and 
DM, arises in the solution of the strong-CP prob¬ 
lem: as in the PQ model, in eq. o 9 is turned into 
a dynamical variable, the phase of $. Here, however, 
the modulus of increases in time, approaching mp 
by our cosmic epoch, when it is DE; meanwhile, 9 
is driven to approach zero, still performing harmonic 
oscillations which are axion DM. The critical time for 
the onset of coherent axion oscillations is the eve of 
the quark-hadron transition, because of the rapid in¬ 
crease of the axion mass m(T, </)). In our dual axion 
scheme, the constant Fpg scale of the PQ model is 
replaced by the slowly varying field (j). Instead of the 
scale FpQ, data fix the scale A, in the SUGRA poten¬ 
tial. This unique setting provides DM and DE in fair 
proportions. 

We therefore simultaneously solve the strong CP 
problem and yield DM and DE in fair proportions by 
setting a single parameter. 
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